Fractal analysis on a closed classical hard-wall billiard using a simplified box-counting 

algorithm 



(N ■ 
O 

o : 

(N ■ 
3 ■ 

<: 
t> : 

(N ■ 



Q 

U 
d 



> 

m 
o 

oo 
o 

O 



x 



Suhan Reef] 

Department of Industrial Information, Kongju National University, 
Yesan-Up, Yesan-Gun, Chungnam, 340-802, South Korea 
(Dated: February 4, 2008) 

We perform fractal analysis on a closed classical hard-wall billiard, the circular billiard with a 
straight cut, assuming there are two openings on the boundary. We use a two-dimensional set of 
initial conditions that produce all possible trajectories of a particle injected from one opening, and 
numerically compute the fractal dimension of singular points of a function that maps an initial 
condition to the number of collisions with the wall before the exit. We introduce a simplified box- 
counting algorithm, which uses points from a rectangular grid inside the two-dimensional set of 
the initial conditions, to simplify the calculation, and observe the classical chaotic properties while 
varying the parameters of the billiard. 
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Chaotic systems have recently attracted many re- 
searchers from various fields, partly because the fast de- 
velopment of computer hardwares has enabled us to solve 
equations almost insolvable in the past[|], For physi- 
cists, the two-dimensional (2D) billiard system has been a 
popular subject for studying the dynamics of the Hamil- 
tonian chaotic systems. Classically, the dynamics of the 
2D billiard system shows three distinct types of behav- 
iors. The system is either integrable (regular behavior) or 
non-integrable (either soft chaos, characterized by mixed 
phase spaces that have both regular and chaotic regions, 
or hard chaos, characterized by ergodicity and mixing) |jj . 

Here we use a closed 2D hard-wall billiard with two 
openings on the boundary, and perform fractal analysis 
using the classical dynamics while varying the shape of 
the bliiard and the size of the openings. To calculate the 
fractal dimension, we use a set of initial conditions that 
will produce trajectories of a particle injected from one 
opening, and calculate certain values for each trajectory, 
such as the exit opening (an opening from which the par- 
ticle exits), the number of collisions with the wall before 
the exit, the dwell time, and so on. When the system 
is chaotic, singular points of this kind of response func- 
tions form a fractal set[f|, |[ ||. In this Letter, we will 
calculate the fractal dimension of these singular points of 
the function that maps an initial condition to the num- 
ber of collisions, by introducing a simplified box-counting 
algorithm. 

The billiard we have chosen is the circular billiard with 
a straight cut, or the "cut-circle" billiard (see Fig. [I]). 
There are five parameters: (1) the width W measured in 
the direction perpendicular to the cut, (2) the radius R, 
(3) the angular width A of the openings, which represents 
the opening size, (4) the orientation angle f2 of the cut 
relative to the first opening, and (5) the position of the 
second opening relative to the first opening as measured 
by the angle 7. We scale the width W by R, so w = W/ R, 



thereby reducing the number of independent parameters 
to four: w, A, fi, and 7. For all subsequent discussions, 
we set f2 = 135° and 7 = 270°. For the closed cut-circle 
billiard (when A = 0), it has been proved in Ref. [7] that 
the phase space is mixed (soft chaos) when < w < 1, 
and that the phase space is fully chaotic (hard chaos) 
when I < w < 2. And it is integrable when w = 1 and 2. 
Consequently this billiard has a single shape parameter w 
that makes the billiard exhibit all three types of chaotic 
behaviors. 

In Ref. ||, a one-dimensional (ID) set of initial condi- 
tions of the cut-circle billiard was used to calculate frac- 
tal dimensions while varying the opening size A and the 
shape parameter w [see Fig. |l](b)] . The particle is in- 
jected with an incident angle (j> (— ir/2 < <fi < tt/2). In 
Ref. authors used a ID set of initial conditions for 
the stadium, and calculated the recurrence time for each 
trajectory to obtain the fractal dimension. These ID 
sets, however, are subsets of the set of all significant ini- 
tial conditions that represent all possible trajectories; in 
our billiard, a 2D subspace in the four-dimensional (4D) 
phase space is a minimal set that will produce all possible 
trajectories originating from one opening [see Fig. |l](c)j. 
(A 2D set of initial conditions were also used in Ref. [3] for 
the Sinai billiard.) There are two reasons that make this 
reduction possible. First, the value of the energy does 
not change the trajectory of the particle. This is due to 
the flat potential inside the billiard with infinitely hard 
boundaries. Second, the distance from the center of the 
billiard to the initial location of the particle can be fixed 
at R. 

With these two kinds of sets, we numerically find a 
function that maps an initial condition to the number of 
collisions before the exit. (In numerical computations, 
all possible initial conditions cannot be used; only a fi- 
nite and discrete subset can be considered.) In Fig. ||, 
we show how the number of collisions changes for these 
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FIG. 1: (a) Geometry of the cut-circle billiard. There are 
two openings, I and II, and a straight cut. The size of the 
cut is given by the width W (here we use w = W/R as the 
relative size of the cut), and the size of two openings is A. 
The positions of the cut and the opening II, Q and 7, are 
measured from the position of the opening I. (b) When a 
particle is injected from the central point of the opening I 
with an angle <f), the set of possible initial conditions is one 
dimensional (—n/2 < <f> < 7r/2). (c) To fully represent all 
significant initial conditions, a particle is injected from a point 
on the circular boundary of the opening I with the angular 
momentum Pgi n - This set, whose element is represented by a 
pair of values, (9i n , Pein), is a closed two-dimensional subspace 
of the four-dimensional phase space. 



sets of initial conditions when w = 1.5 and A = 60°. 
In Fig. ||(a), the graph of 4> vs the number of collisions 
is shown when the ID subset of initial conditions as in 
Fig. 0(b) is used (see Ref. || for more results). In Fig. 
H(b), the graph shows the full 2D set of initial conditions, 
(#i„, P$i n ), bounded by two equations, 
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FIG. 2: Calculations of the number of collisions when 
n = 135°, 7 = 270°, w = 1.5, and A = 60°. (a) Graph 
of the incident angle (f> vs the number of collisions Nc when 
the particle is injected as in Fig. [j](b). (b) When the parti- 
cle is injected as in Fig. 0(c), all significant initial conditions 
lie inside the closed two-dimensional space bounded by the 
bold line in the graph. In the graph, only singular points are 
shown. The dotted line represents the one-dimensional sub- 
set of initial conditions used in (a) . Five geometrical channels 
with low Nc's are also shown. 



which are represented by the bold line. Instead of the 
number of collisions, the graph shows singular points that 
constitute singular boundaries at which the number of 
collisions changes. We can observe that singular points 
have infinitely fine structures. The number of collisions 
are shown only for five geometrical cftarmefe jH], ^) in 
this example. The dotted line is the ID subset used in 
Fig. |(a). 

Now we calculate the fractal dimension df of a set of 
singular points in this 2D set using a simplified box- 
counting algorithm. As a generalization of the box- 
counting algorithm used in Ref. ||for the ID set, we use 
a uniform 2D rectangular grid. Then there are uniform 
rectangular boxes, and each of them is represented by 
four points on the grid. For all of these uniformly dis- 
tributed grid points inside the 2D set, we numerically 
find the number of collisions before the exit, and for each 
box, we compare the number of collisions at four corners. 
With an assumption that a box does not contain a singu- 
lar point when four values are all equal, we can count the 
number of boxes containing singularities out of all boxes 
inside the 2D set. The number of all boxes inside the set 
is Np, and the number of boxes that contain any singular 
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FIG. 3: (a) When w — 1.5, the number of boxes containing 
the singular points, Ns, is counted as the number of boxes, 
Np, increases, for five different A values. Then df can be 
found from the slopes of these log-log graphs, (b) Graphs of 
df vs the opening size A for w — 0.75, 1.04, and 1.5. Here we 
compare the values of df for the 2D set (o) and for the ID 
set (*). 



point, based on the above assumption, is N$- Then the 
fractal dimension df is denned by 



lim 

N P 



logio N s 
°° log 10 Np 1 ' 2 ' 



(3) 



which is the slope of the graph of 0.51og 10 /Yp vs 
log 10 Ns- In numerical calculations, we find N$ for sev- 
eral Np values, and then use the ordinary least-square 
fit to find the slope using the points in the graph. 

In Fig. ||(a), the log-log graph is shown for several A 
values when w — 1.5. (For all calculations, Np up to 10 8 
will be used.) Unlike the calculations using the ID set, 
df will be in the range of 1 < df < 2. We can define df, 
which better represents the fractality, 



d f = 



2-df (for the 2D set) 
1 - df (for the ID set) 



(4) 
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FIG. 4: Graphs of w vs d/ for seven different opening sizes: 
A = 1° (o), A = 5° (□), A = 10° (O), A = 20° (A), A = 30° 
(V), A = 40° (+), and A = 50° (*), for the 2D set. When 
A is small enough (A < 20° in the graph), the behavior of 
the graph is clearly distinct for two regions: hard chaos (1 < 
w < 2) and soft chaos (0 < w < 1). But for larger A, the 
distinction between two regions starts to disappear. 



In Fig. |(b), df's both for the 2D set and for the ID set 
are compared for three different w values while varying 
A. Like the curves for the ID set, not all curves mono- 
tonically increase as A increases, even though it is not 
as significant. The reason for this phenomenon is that 
there are more possible trajectories when the opening 
gets bigger, and it also explains why curves don't match 
very well for large A. We also observe that the cases 
with hard chaos (w — 1.04, and 1.5) do not behave pre- 
dictably for large A, which will be more clearly seen in 
the following graph. When the original shape is changed 
sufficiently by the presence of the openings, the ergodic- 
ity of the billiard is no longer an important factor; the 
relative locations of the openings start to have more effect 
on the fractal dimension as in the cases of soft chaos. 

In Fig. |[ we calculate df as the width w varies from 
0.5 to 2 with a step size of 0.02 for seven opening sizes 
(A = 1°, A = 5°, A = 10°, A = 20°, A = 30°, A = 
40°, and A = 50°), using the 2D set. We can compare 
graphs in two different regions (found in the closed cut- 
circle billiard): < w < 1 (soft chaos) and 1 < w < 
2 (hard chaos). We observed similar behaviors seen in 
results for the ID set in Ref. pi However df does not 
reach one (df = 1 represents non-chaos), even when the 
billiard is integrable (w = 1 or 2), and this is due to 
the fine structure near trajectories bouncing very close 
to the circular boundary. When the opening size A is 
not big (see graphs for A = 5°, A = 10°, and A = 20°), 
we observe that the behavior in two regions is clearly 
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distinct. In the region < w < 1, there are fluctuations, 
which comes from the mixed phase space structures of the 
billiard, and in the region 1 < w < 2, graphs are smooth 
because the phase spaces of the billiard have no structure 
due to ergodicity. On the other hand, when the opening 
size gets bigger (see graphs for A = 30°, A = 40°, and 
A = 50°), the distinction between two regions, observed 
in cases with smaller openings, disappears. There are 
fluctuations in both regions. 

In this Letter, we calculated the fractal dimension for 
the cut-circle billiard with two openings by introducing 
the simplified box-counting algorithm for the 2D set of 
initial conditions. The simplified box-counting algorithm 
introduced here is possible because the function used to 
find singularities is constant except singular points (i. e., 
the function is constant with different values in all regions 
surrounded by singular boundaries). With this method, 
some boxes containing singular points will be missed, but 
the number of missed boxes is negligible when we calcu- 
late the fractal dimension with big enough Np values. 
Results for the 2D set were close to those for the ID set 
only when A is small enough; hence using the ID set 
in fractal analysis does not fully represent the chaoticity 
of the billiard for most cases. In conclusion, this kind 
of fractal analyses gives us one of the fundamental mea- 
sures of the chaoticity of classical billiard systems, and 
finding the relations with the quantum and semiclassical 
dynamics of the same kind of billiards will be one of the 



interesting future works. 
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